Universal Control of Nuclear Spins Via Anisotropic Hyperfine Interactions 
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We show that nuclear spin subsystems can be completely controlled via microwave irradiation of 
resolved anisotropic hyperfine interactions with a nearby electron spin. Such indirect addressing of 
the nuclear spins via coupling to an electron allows us to create nuclear spin gates whose operational 
time is significantly faster than conventional direct addressing methods. We experimentally demon- 
strate the feasibility of this method on a solid-state ensemble system consisting of one electron and 
one nuclear spin. 



Coherent control of quantum systems promises opti- 
mal computation jl] , secure communication [2] , and new 
insight into the fundamental physics of many-body prob- 
lems [3]. Solid-state proposals [H |5J IHl El H] for such 
quantum information processors employ isolated spin de- 
grees of freedom which provide Hilbert spaces with long 
coherence times. Here we show how to exploit a local, 
isolated electron spin to coherently control nuclear spins. 
Moreover, we suggest that this approach provides a fast 
and reliable means of controlling nuclear spins and en- 
ables the electron spins of such solid-state systems to be 
used for state preparation and readout j9] of nuclear spin 
states, and additionally as a spin actuator for mediating 
nuclear-nuclear spin gates. 

Model System. The spin Hamiltonian of a single 
local electron spin with angular momentum, = j and 
N nuclear spins, each with angular momentum Ik = \, 
in the presence of a magnetic field B is [TD] : 
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Here /3e is the Bohr magneton, 7^ is the gyromagnetic 
ratio; S and I'' are the spin-| operators. The second- 
rank tensors g, A'^, S, and D represent the electron g- 
factor, the hyperfine interaction, the chemical shift, and 
the nuclear dipole-dipole interaction respectively. 

In the regime where the static magnetic field B ^ Bqz 
provides a good quantization axis for the electron spin, 
the Hamiltonian can be simplified by dropping the non- 
secular terms which corresponds to keeping only electron 
interactions involving Sz- The quantization axis of any 
nuclear spin depends on the magnitudes of the hyperfine 
interaction and the main magnetic field, as well as their 
relative orientations. When these two fields are compa- 
rable in magnitude |37) J^q can be approximated by: 
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FIG. 1: Energy level diagram for Eq. ^ with N=l. The 
electron spin state is in an eigenstate of purely the Zeeman 
interaction, while the nuclear spin state is not an eigenfunc- 
tion of the Zeeman interaction alone due to the anisotropic 
hyperfine interaction. Because (ao|/3i) 7^ and {ao|/3o) 7^ 
the electron spin operator (Sx) has finite probabilities between 
all levels (dashed arrows). This allows for universal control of 
the entire spin system. The filled and unfilled circles represent 
the relative spin state populations of the ensemble at thermal 
equilibrium. In our experimental setup the energy differences 
are Loi2/2n = 7.8 MHz, 1.^34/271 = 40 MHz, a;i4/27r = 12.005 
GHz, tJ23/27r = 11.954GHz 
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The nuclear dipole-dipole interaction is neglected as it is 
typically 10^ times weaker than the hyperfine terms. 

As described in Figure [l] (N=I), the nuclear spin is 
quantized in an effective field that is the vector sum of 
the Zeeman and hyperfine interaction. The nuclear spin 
eigenstates are a mixture of the nuclear Zeeman eigen- 
states: |ao) 
to lao). 
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An N-bit string of quantum information can be stored 
in this system in either the spin-up or spin-down electron 
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manifolds: 

\ji 32 ■ .-Jn) = IT) ® WD ® Wl) • • • |af„) (3) 
= |i)®|/3j,)®|/3^)0---®|/3f„) (4) 

where ju G {0,1}. Note that by storing information in ei- 
ther the It) or ID electron spin manifolds there is no spin 
superposition of the electron wavefunction and that the 
electron spin remains separable from the nuclear spins. 

Universal Control of 2'^"'"^-dimensional Hilbert 
Space. Having identified a system for storing quan- 
tum information wc show the condition under which such 
a system allows for arbitrary quantum operations. A 
necessary and sufficient condition for complete control- 
lablity, or universality, requires the nested commutators 
of the natural Hamiltonian Eq. ^ and each of the 
control Hamiltonians are a closed Lie group SU(2^+^) 
[38j I12j . An equivalent diagrammatic representation 
relies on graph connectivity for assessing the controlla- 
bility of quantum systems represented as Lie algebras 

[131 [n]. 

Our controls are time-dependent microwave fields os- 
cillating at the electron spin resonance frequency {'^ 
PedzzBo) and parametrized by three values: Bf, the am- 
plitude of an oscillating magnetic field (_L Bqz), the 
frequency of oscillation, and the phase of the oscilla- 
tion: 

/3eg,,Bf(t)(cos(r!i + 0(t)) + sin + 4) 

(5) 

Barring degenerate eigenvalues Ej of or degenerate 
transition frequencies {hujjk = Ej — Ek), the matrix rep- 
resentation of is strongly regular. As shown in |13| . 
this and the complete connectivity of the graph gener- 
ated by the matrix elements of J^c guarantee universal- 
ity. For the le-N nuclear spin system, distinct gyromag- 
netic ratios and hyperfine couplings for each nuclear spin 
guarantee the non-degeneracy of the eigenstates. The hy- 
perfine couplings and the Zccman frequencies must also 
be chosen such that LUj^/LOju.' =/= 1. Lastly, the anisotropy 
of each hyperfine interaction assures the complete con- 
nectivity of the graph given the form of J>^o and J^c- 
Figure [2] shows several diagrammatic examples. For uni- 
versal control over a set of nuclear spins with resolved 
anisotropic hyperfine coupling to one electron spin, it is 
sufficient to apply only an amplitude modulated wave- 
form to any electron spin transition at a fixed f2. 

We can provide some insight into how modulation of 
the electron spin state through shaped microwave fields 
provides control over the nuclear spins. Flipping the elec- 
tron spin changes the quantization axis of the nuclear 
spin. Since the these two quantization axes are separated 
by a large angle, a sequence of spin evolutions under these 
two non-commuting axes permits arbitrary nuclear rora- 



tions. If we consider collective motions of the N nuclear 
spins relative to the two electron spin states, this gener- 
ates the complete algebra in SU(2^"'"^). Given the com- 
plexity of the full dynamics of these 2^ vectors, we can 
use optimal control methods developed for and applied to 
liquid state nuclear magnetic resonance (NMR) to engi- 
neer arbitrary unitaries [13 [121 [H] • We limit the control 
fields of the electron-nuclear system to only the electron 
spin flip transitions and achieve quantum gates whose 
operation times are faster than if we relied upon nuclear 
spin nutation. The idea of using pulsed ESR to quickly 
generate nuclear spin coherence and echoes was flrst used 
for spectroscopy [HI [19]; here, we extend this idea to 
not only generate nuclear coherence, but also suppress 
all other closed-system dynamics to generate a desired 
quantum gate. For organic crystal systems of a few nu- 
clear spins the values of the hyperfine interaction are tens 
of MHz. For pulsed ESR systems the Rabi frequencies 
for an electron spin are also tens of MHz [T0 | [20 | |2T | |22] . 
while the nuclear Rabi frequency is at most 1 MHz |23j . 
The nuclear- nuclear dipolar coupling is tens of kHz. 
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FIG. 2: The connectivity of energy levels of le-Nn systems 
represented as 2^"'"'^ node graphs. An edge is drawn between 
two nodes if the control Hamiltonian operator has a non-zero 
matrix element for the eigenstates represented by the nodes 
{{k\Sx\j) 7^ 0). When the hyperfine interaction between any 
one nuclear spin and the electron is purely isotropic univer- 
sality is not achieved as in (a). 



The le-ln system. We demonstrate the utility of 
this control scheme by exploring Ramsey fringes |24j 
and Hahn echoes [2S] in a le-ln system — a single 
crystal of x-ray irradiated malonic acid [3BJ [271 
The parametrized Hamiltonian is J%/2tt = v^Sz — 
yJz + AzxSzL + AzzSjz where i^s = 11.885 GHz, 
= 18.1 MHz, A^^j, « 14.2 MHz, and A^^ « -42.7 
MHz. The control Hamiltonian Eq. ^ parameters are 
max(/3eg^^Sf(t))/27r = 7 MHz and n/2Tr = 11.909 GHz. 
No arbitrary phase controls cjj^t) were used. All experi- 
ments were performed on a home-built pulsed ESR spec- 
trometer with a 12W solid-state microwave amplifier ca- 
pable of amplitude modulations of 250 Ms/s using an ar- 
bitrary waveform generator (AWG) and double-balanced 



3 



mixer. A home-built low-temperature probehead cools 
the sample to 77K and contains a loop gap resonator 
with a volume of ~ 55mm'^ and a Q of 250. 

Experimental Results. To implement an arbitrary 
unitary propagator we use the GRadient Ascent Pulse 
Engineering (GRAPE) [TS] algorithm for finding the con- 
trol field J^Q. Constraints on the the modulation se- 
quence, such as maximum nutation rate and pulse band- 
width, were chosen in accordance with our hardware lim- 
itations such as finite power amplifiers, modest AWGs, 
and finite bandwidth components. The simulated gate 
fidelity jTlj of useful gates is at least 0.98. 

The equilibrium state of the ensemble system, 
/Othcrmai = —Sz, has no net nuclear spin polarization, so 
we first transfer the available electron spin polarization 
to the nuclear spins. This is achieved by selectively in- 
verting the levels |2) and |4) or |1) and |3) (see Figure [T]). 
We created coherence between nuclear eigenstates with 
an engineered nuclear 7r/2 pulse selective for only one 
of the electron manifolds, Ui2{'k/2) [35]. The Ramsey 
fringe experiment measures the phase evolution under 
Mq. We halt evolution by again applying Ui2{tt/2) and 
then transferring the polarization back to the electron 
spin. By monitoring the relative amplitude of the elec- 
tron spin echo at different times, r, we indirectly observe 
the nuclear spin dynamics. If we introduce a refocusing 
pulse for both nuclear spin states, Ur = Ui2{n) © [734(77), 
after a time r/2, the coherent phase oscillations are re- 
focused. 



Figure [4] shows the coherent oscillations between nu- 
clear coherence on levels |1) & |2). We implement the 
net unitaries, Upc = Ui2{^)U24{tt) and as a single 
modulation sequence with total time (T) of 800ns and 
simulated fidelities (F) 0.99 and 0.98 respectively. The 
nuclear tt pulse, Ur, was implemented in a time period 
(F=0.98, T=520ns) much shorter than would have been 
possible by addressing the nuclear spin transitions di- 
rectly. Again, it is key to realize that the nuclear pulse is 
achieved through modulation of the hyperfine interaction 
and is applied at the electron spin resonance frequency. 
Simulations of the modulation sequences using our model 
le-ln Hamiltonian show agreement of the observed oscil- 
lation. 
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FIG. 3: Color online. Schematic pulse sequence for mea- 
suring Ramsey fringes and Hahn echoes. Up creates a non- 
equilibirum population difference between levels 1&2 and Uc 
creates a coherence between nuclear spin in the S=-l/2 man- 
ifold. During r, this coherence evolves under J^, acquiring 
an observable phase. The coherence is transformed back to 
a population difference between nuclear spin levels and then 
to electron spin levels. A pair of short, unmodulated pulses 
are used to detect an electron spin echo whose height is pro- 
portional to the resultant electron spin population. With a 
refocusing pulse (?7r) the acquired phase is unraveled, leav- 
ing no modulation of the echo signal. The waveform used to 
implement UpUc is shown inset. Note that all pulses are ap- 
plied resonant with the 1-4 transition and induces transitions 
between 1-4, 2-4, 1-3, and 2-3 due to selection rules. 



FIG. 4: Color online. Measurements of the electron spin echo 
as a function of r between coherence transfer indirectly reveal 
the nuclear precession rate. Numerical simulations of the ex- 
periment show agreement with the observed signal (signal-to- 
noise ratio « 4.9). The Ramsey fringe experiment (x) reveals 
a clear precession of the coherence between the |1) and |2) 
states at roughly 8 MHz . For the echo data (□), small pulse 
errors and incomplete phase cycling result in a small system- 
atic oscillation at half the precession frequency; ideally the 
echo should be independent of r. 



Conclusions and Outlook. Using shaped pulse 
techniques for engineering quantum gates, we have de- 
scribed a method to control nuclear spins coupled via the 
hyperfine interaction to a nearby localized electron spin 
impurity using only electron spin transitions. Further- 
more we have experimentally demonstrated the viability 
of these ideas using a testbed system of one electron spin 
(S=l/2) and one nuclear spin (1=1/2). 

Using the same experimental setup this method can, 
in principle, extend to any number of nuclear spins 
that have resolvable anisotropic hyperfine interactions. 
Specifically, we have found a le-2n CNOT gate, (see 
supporting materials for details) between the carbon and 
proton nuclear spins for malonic acid isotopically labeled 
with a ^■^C at the methylene position This gate 

can be performed in only 2/is. If only nuclear nuta- 
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tion frequencies and nuclear dipole-dipole interactions 
were used, such a CNOT would take much longer due 
to the relative strength of D to A. Modulations of the 
anisotropic hyperfine interactions can generate quantum 
gates between nuclear spins without changing the state 
of or entangling the electron spin. Previously studied 
solid-state spin systems with multiple resolved nuclear 
spins O [30l im mi |33] can thus be controlled by only 
modulating the electron spin states. Moreover, it should 
be possible to augment the controllable Hilbert space by 
coupling localized electron spin states. Such coupling 
may be electrical, with systems like phosphorous-doped 
silicon |34j or optical in systems akin to nitrogen-vacany 
defects in diamond [35] . As the ultimate utility of these 
systems depends on their decoherence, characterization 
of these mechanisms is key. Precise unitary engineering 
can also be used to measure relaxation processes or per- 
form quantum process tomography in hyperfine coupled 
systems in a more controlled manner than has been re- 
ported [2D]. 
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Local 

Zeeman + Hyperfine = Nuclear 
Field 



FIG. 5: Contributions to nuclear field. The nuclear Zeeman field adds to the local hyperfine field yielding the quantization axis 
for the nuclear spin. Because the electron spin can be aligned or anti-aligned with the main magnetic field, the local hyperfine 
field can be either positive or negative. The nuclear spin aligns (or anti-aligns) with the local magnetic field. The presence 
of an anisotropic term hyperfine interaction comparable to the nuclear Zeeman interaction prevents the nuclear spins aligning 
parallel or antiparallel to the electron spins. Note that the quantization axis of the nuclear spin differs by a finite angle in these 
two cases. 
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Electron-Nuclear SWAP gate, Malonic Acid 
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FIG. 6: Electron-Nuclear SWAP gate. For the same field and crystal orientation used in the above experiment, we find a 
modulation sequence (a) to perform a SWAP gate between the electron and the nuclear spin. The sequence consists 420 4ns 
intervals each specifying a pulse amplitude from -1 to 1, corresponding to ±ljI the electron spin Rabi frequency. This gate 
has a simulated fidelity of 0.999 and has been experimentally implemented. The frequency domain (b) shows the spectral 
components of the sequence. Note that nearly half of the power is contained outside the ±25MHz bandwidth centered about 
0. If we take into account the Q of the resonator, thus attenuating these higher frequency components, the simulated fidelity 
degrades to 0.964. 
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FIG. 7: Nuclear-Nuclear CNOT gate. If the carbon nucleus closest to the paramagnetic defect of malonic acid is a C isotope 
(1=1/2), we have a le-2n system. Using hyperfine tensor values reported in the literature, we find a modulation sequence to 
implement a CNOT gate between two nuclear spins: i7ie-2n = 1 (8) ChNOTc- At the same crystal orientation used above, 
we find a sequence modulating both the amplitude and the phase of the control Hamiltonian with 1000 2ns intervals and a 
ma:ximum Rabi frequency of 15MHz. The time for performing the gate using nuclear-nuclear dipole couplings (tens of kHz) in 
the absence of an electron spin would be much longer. In our scheme the nuclear-nuclear dipole couplings enter as a source of 
decoherence and can be included in the model to generate gates robust to this interaction. 
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FIG. 8: Trajectories of spin systems on projected Bloch Spheres. After a polarization transfer from the equihbrium electron 
spin state to the nuclear spin, we apply a 90° rotation between the |1) and |2) manifold. The trajectory of the spin vector can 
be traced by taking an equivalent Bloch sphere between any two levels, where the usual operators {ax,(Jy, ctz} are analogous to 
<^^* = b)(^l + l'!^)OI) '^i'^ ~ and a^J^ — \ — For the 12-pair, the state begins along the z axis (green) 

and then rotates down to the equatorial plane (red). The 14-pair reduces it's z component by a factor of two, as half of the 
state goes into coherence; the 34-pair begins and ends along the z-axis, as the rotation is in a subspace orthogonal to 34. 
Although the trajectories the spin takes can be very complicated, only the starting and end points are significant. 



